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DESCRIPTION 



SAMPLING RATE CONVERTER AND METHOD OF SAME AND AUDIO 

APPARATUS 

5 

TECHNICAL FIELD 
[0001] The present invention relates to a sampling 
rate converter which can be applied to for example 
conversion of a sampling frequency of audio and 

10 conversion of resolution for magnification or reduction 

of a frame of an image and a method of the same and audio 
apparatus . 

BACKGROUND ART 
[0002] In digital signal processing of audio and 

15 images, filter processing is frequently used. As the 

filter used for the filter processing, due to the feature 
that a linear phase is obtained with a finite number of 
taps, a linear phase FIR (Finite Impulse Response) filter 
is frequently utilized. 

20 [0003] FIG. 1 is a diagram showing the configuration 
of a transversal type circuit of the linear phase FIR 
filter. 

This linear phase FIR filter 1 has, as shown in FIG. 
1, (n-1) number of delay units 2-1 to 2-n-l cascade 
25 connected with respect to an input terminal TIN and 

configuring a shift register, n number of multipliers 3-1 
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to 3-n for multiplying filter coefficients h(0) to h(n-l) 
with respect to the signal input to the input terminal 
TIN and output signals of the delay units 2-1 to 2-n-l, 
and an adder 4 for adding the output signals of the n 
5 number multipliers 3-1 to 3-n and outputting the result 
to an output terminal TOUT. 

[0004] As a representative design method of such a 
linear phase FIR filter, for example a Remez Exchange 
algorithm applied to a linear phase FIR filter by Parks, 
10 T.W. and McClellan, J.H. et al . is known (refer to for 
example Non-Patent Document 1) . 

[0005] The Remez Exchange algorithm is an algorithm 
for approximation so that a weighted approximation error 
exhibits an equal ripple shape with respect to a desired 

15 amplitude characteristic. 

[0006] As applications of filter processing using a 
linear phase FIR filter, there are conversion of 
resolution of an image utilizing sampling rate conversion 
and conversion of the sampling frequency of audio. 

20 For example, in conversion of resolution, use is 

made of a multi-rate filter using an interpolator, a 
decimeter, and a linear phase FIR filter as element 
technologies (refer to for example Non-Patent Document 2) . 
[0007] In a multi-rate filter, generally the linear 

25 phase FIR filter is used poly-phase decomposed matching 
with the interpolator. Both of the interpolator and the 
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decimeter are cyclical time-invariant systems and have 
characteristic features different from that of a time- 
i nvar i an t sy s tern . 

Due to the cyclical time invariability of the 
5 interpolator, distortion on a lattice called as 

"chessboard distortion" occurs in the conversion of 
resolution of the image. 

[0008] Therefore, Harada and Takaie considered 
conditions for avoiding the chessboard distortion from a 
10 zero point arrangement of the filter (refer to Non-Patent 
Document 3) . 

[0009] A transmission function H(z) of the multi-rate 
filter not accompanied by chessboard distortion is found 
by multiplying a transmission function K(z) of the linear 
15 phase FIR filter (hereinafter referred to as an 

equalizer) designed by some sort of method by a zero 
point transmission function Z(z) in order to avoid the 
chessboard distortion later. 
[0010] (Equation 1) 
20 H(z)=Z(z) K(z) (1) 

[0011] (Equation 2) 

Z (z)=l + z" 1 +Z" 2 + +Z" (u " 1) " (2) 

[0012] Here, a previously fixed linear phase FIR 
filter like the zero point transmission function Z(z) for 
25 avoiding the chessboard distortion will be referred to as 
a pre-filter. 
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[0013] FIGS. 2A to 2C show an example of the frequency 
response of a multi-rate filter avoiding chessboard 
distortion by multiplying the equalizer designed by the 
Remez Exchange algorithm by a pre-filter and the weighted 
5 approximation error. 

Non-Patent Document 1: Parks, T.W. and McClellan, 
J.H.: "Chebyshev Approximation for Nonrecursive Digital 
Filters with Linear Phase", IEEE Trans. Circuit Theory, 
CT-19, 2, pp. 189-194, 1972, and Rabiner, L.R. , McClellan, 

10 J.H. and Parks, T.W. : "FIR Digital Filter Design 

Techniques Using Weighted Chebyshev Approximation", Proa. 
TEEE, Vol. 63, April, pp. 595-610, 1975; 

Non-Patent Document 2: Takaie, Hitoshi, Multi-rate 
Signal Processing, Shokodo, 1997; 

15 Non-Patent Document 3: Harada, Yasuhiro and Takaie, 

Hitoshi: Multi-rate Filter not Accompanied by Chessboard. 
Distortion and Zero Point Arrangement of Same, Shingaku 
Giho CAS96-78, pp. 1-6, 1997-01 

DISCLOSURE OF THE INVENTION 

20 PROBLEM TO BE SOLVED BY THE INVENTION 

[0014] However, methods for avoiding chessboard 
distortion by the above-described methods has the 
following disadvantages. 

Namely, in a multi-rate filter having a 

25 transmission function H(z) designed by the conventional 
method, as shown in FIG. 2C, the equal ripple of the 
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weighted approximation error designed by the Remez 
Exchange algorithm collapses. 

Further, in the multi-rate filter designed by the 
conventional method, as shown in FIG . 2B, the pass band 
5 gain is not a constant value, but the right end 
attenuates * 

[0015] When conversion of resolution is carried out by 
using such a filter, the contours of an image become 
vague and influences the image quality. In the same way, 

10 when conversion of the sampling frequency is carried out 
by using such a filter, a high frequency component 
attenuates and high precision audio reproduction cannot 
be carried out. 

Attenuation at the pass band cannot be avoided even 

15 when the filter coefficient is increased. 

[0016] Further, there is a case where a DC gain must 
be set at 1 when the frequency w is 0 in the design 
specifications of a linear phase FIR filter. 

However, in a Remez Exchange algorithm, as shown in 

20 FIGS. 3A to 3C, an amplitude characteristic that passes 
any designated frequency point cannot be obtained. 
[0017] Further, FIG. 4 shows an example of the 
configuration of a sampling rate conversion unit based on 
the theory of conventional multi-rate signal processing. 

25 A sampling rate conversion unit 10 of FIG. 4 has an 

up sampler 11, FIR filter 12, and down sampler 13. 
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Further, in FIG. 4, "U" and "D" are positive 
integers primary with respect to each other, and H(z) 
indicates the transmission function of the FIR filter. 
Further, an upward arrow indicates the up sampler for 
5 inserting (U-l) zero points between signals, and a 

downward arrow indicates the down sampler for thinning 
out signals at an interval of D signals. 

[0018] The sampling rate conversion unit 10 of FIG. 4 
raises the sampling rate of the signal U-fold by the up 
10 sampler 11, performs band restriction by using the FIR 

filter 12, and finally performs an operation lowering the 
sampling rate to 1/D by the down sampler 13. 

By this, the sampling rate can be converted to U/D 
times. The cutoff frequency of this FIR filter becomes as 
15 follows. 

[0019] (Equation 3) 

co C= 7i /maxU , D ( 3 ) 

[0020] For example, in order to accomplish sampling 
rate conversion from 44.1 kHz or up to 44.1 kHz based on 
20 the theory of the conventional multi-rate signal 

processing as shown in FIG. 4, a filter in which the 
cutoff frequency of the FIR filter 2 of the low pass 
filter is extremely strict must be prepared. 

In general, as shown below, such a filter needs a 
25 large number of taps and is extremely difficult to 
realize. 

6 
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[0021] • From 44.1 kHz to 48 kHz, 
Up=160 / Down=147, and 
Cutoff frequency 7t/160. 

• From 44.1 kHz to 32 kHz, 

Up=320, Down=441 , and 

Cutoff frequency 7t/441. 
■ From 44.1 kHz to 24 kHz, 
Up=80, Down=147, and 
Cutoff frequency 7c/147. 

• From 48 kHz to 44.1 kHz, 

Up=147, Down=160, and 

Cutoff frequency 7r /160 . 

• From 32 kHz to 44.1 kHz, 

Up=4 4 1 , Down=32 0 , and 
Cutoff frequency 7t/441. 

• From 24 kHz to 44.1 kHz, 

Up=147, Down=80, and 

Cutoff frequency 7c/147. 
[0022] The present invention was made in consideration 
with such a circumstance and has an object thereof to 
provide a sampling rate converter in which the equal 
ripple of the weighted approximation error does not 
collapse, the gain of the pass band can be held at 
substantially a constant value, an amplitude 
characteristic that passes any frequency can be obtained, 
and high precision conversion can be accomplished and a 
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method of the same and an audio apparatus . 

MEANS FOR SOLVING THE SUBJECT 
[0023] To attain the above object, a sampling rate 
converter according to a first aspect of the present 
5 invention has an up sampler for inserting U-l zero points 
between sample signals and raising a sampling frequency 
U-fold, a convolution processing unit including an FIR 
filter and performing predetermined convolution 
processing with respect to an output signal of the up 

10 sampler, and a linear interpolation block for selecting 
two points of samples with respect to the results of 
processing of the convolution processing unit and finding 
a value at a required position from the linear 
interpolation, wherein the FIR filter of the convolution 

15 processing unit is an FIR filter where an impulse 
response is expressed by a finite time length, the 
impulse response becomes a filter coefficient, and a 
transmission function H(z) is associated with a 
transmission function Z(z) of a pre-filter, and the 

20 filter coefficient is set by performing weighted 

approximation with respect to a desired characteristic in 
relation to a frequency to be passed and/or a frequency 
response of the pre-filter. 

[0024] A sampling rate converter according to a second 
25 aspect of the present invention has a plurality of 

convolution processing units including pre-phase filters 
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obtained by poly-phase decomposing a predetermined FIR 
filter and performing the convolution processing of input 
sample signals and the poly-phase filters decomposed to 
the poly-phases, a plurality of up samplers for inserting 
5 U-l zero points between output signals of corresponding 
convolution processing units and raising the sampling 
frequency U-fold, an adding means for generating a signal 
after adding all signals by adjusting a propagation time 
of output signals of the plurality of up samplers, and a 

10 linear interpolation block for selecting two points of 
samples with respect to the signal by the adding means 
and finding the value at the required position from the 
linear interpolation, wherein the FIR filter is an FIR 
filter where an impulse response is expressed by a finite 

15 time length, the impulse response becomes the filter 
coefficient, and a transmission function H(z) is 
associated with a transmission function Z(z) of a pre- 
filter, and the filter coefficient is set by performing 
weighted approximation with respect to a desired 

20 characteristic in relation to a frequency to be passed 
and/or a frequency response of the pre-filter. 
[0025] A sampling rate converter according to a third 
aspect of the present invention has a convolution 
processing unit including poly-phase filters able to set 

25 different filter coefficients obtained by poly-phase 
decomposing a predetermined FIR filter and performing 
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convolution processing of input sample signals and a 
poly-phase filter having a selected coefficient, a 
selector for selecting two points of samples required for 
an output sample and selecting the coefficient of the 
5 corresponding poly-phase filter, and a linear 

interpolation block for finding the value at the required 
position from linear interpolation, wherein the FIR 
filter is an FIR filter where an impulse response is 
expressed by a finite time length, the impulse response 

10 becomes the filter coefficient, and a transmission 

function H(z) is associated with a transmission function 
Z(z) of a pre-filter, and the filter coefficient is set 
by performing weighted approximation with respect to a 
desired characteristic in reflation to a frequency to be 

15 passed and/or a frequency response of the pre-filter. 

[0026] A sampling rate conversion method according to 
a fourth aspect of the present invention comprising a 
first step of inserting U-l zero points between sample 
signals and raising the sampling frequency U-fold, a 

20 second step of performing predetermined convolution 

processing with respect to a signal multiplied in its 
sampling frequency by U by a convolution processing unit 
including an FIR filter in which an impulse response is 
expressed by a finite time length, an impulse response 

25 becomes the filter coefficient, and a transmission 

function H(z) is associated with a transmission function 
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Z(z) of a pre-filter, and a third step of selecting two 
points of samples with respect to the results of 
processing and finding the value at the required position 
from linear interpolation, wherein the filter coefficient 
5 of the FIR filter is calculated by performing weighted 

approximation with respect to a desired characteristic in 
relation to a frequency to be passed and/or a frequency 
response of the pre-filter. 

[0027] A sampling rate conversion method according to 

10 a fifth aspect of the present invention comprising a 

first step of performing convolution processing of input 
sample signals and poly-phase filters decomposed to poly- 
phases by a plurality of convolution processing units 
including poly-phase filters obtained by poly-phase 

15 decomposing a predetermined FIR filter, a second step of 
inserting U-l zero points between output signals of 
corresponding convolution processing units and raising 
the sampling frequency U-fold, a third step of adjusting 
the propagation time of a plurality of signals having 

20 sampling frequencies raised U-fold and generating a 

signal obtained by adding all signals, and a fourth step 
of selecting two points of samples with respect to the 
signal by the third step and finding the value at the 
required position from the linear interpolation , wherein 

25 the FIR filter is a FIR filter where an impulse response 
is expressed by a finite time length, the impulse 
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response becomes the filter coefficient, and a 
transmission function H(z) is associated with a 
transmission function Z(z) of a pre-filter, and the 
filter coefficient is calculated by performing weighted 
5 approximation with respect to a desired characteristic in 
relation to a frequency to be passed and/or a frequency 
response of the pre-filter. 

[0028] A sampling rate conversion method according to 
a sixth aspect of the present invention comprises a first 

10 step of selecting two points of samples required for an 
output sample and selecting a coefficient of a 
corresponding poly-phase filter and a second step of 
performing convolution processing of input sample signals 
and the poly-phase filter having the selected coefficient 

15 by a convolution processing unit including poly-phase 
filters obtained by poly-phase decomposing a 
predetermined FIR filter and able to set different filter 
coefficients, wherein the FIR filter is a FIR filter 
where an impulse response is expressed by a finite time 

20 length, the impulse response becomes the filter 
coefficient, and a transmission function H(z) is 
associated with a transmission function Z (z) of a pre- 
filter, and the filter coefficient is calculated by 
performing the weighted approximation with respect to a 

25 desired characteristic in relation to a frequency to be 
passed and/or a frequency response of the pre-filter. 
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[0029] A seventh aspect of the present invention is an 
audio apparatus including a sampling rate converter, 
wherein the sampling rate converter has an up sampler for 
inserting U-l zero points between sample signals and 
5 raising a sampling frequency U-fold, a convolution 

processing unit including an FIR filter and performing 
predetermined convolution processing with respect to an 
output signal of the up sampler, and a linear 
interpolation block for selecting two points of samples 

10 with respect to the results of processing of the 

convolution processing unit and finding a value at a 
required position from linear interpolation, wherein the 
FIR filter of the convolution processing unit is an FIR 
filter where an impulse response is expressed by a finite 

15 time length, the impulse response becomes the filter 
coefficient, and a transmission function H(z) is 
associated with a transmission function Z(z) of a pre- 
filter, and the filter coefficient is set by performing 
weighted approximation with respect to a desired 

20 characteristic in relation to a frequency to be passed 
and/or a frequency response of the pre-filter. 
[0030] An eighth aspect of the present invention is an 
audio apparatus including a sampling rate converter, 
wherein the sampling rate converter has a plurality of 

25 convolution processing units including pre-phase filters 
obtained by poly-phase decomposing a predetermined FIR 
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filter and performing convolution processing of input 
sample signals and poly-phase filters decomposed to poly- 
phases, a plurality of up samplers for inserting U-l zero 
points between output signals of corresponding 
5 convolution processing units and raising the sampling 

frequency U-fold, an adding means for generating a signal 
after adding all signals by adjusting a propagation time 
of output signals of the plurality of up samplers, and a 
linear interpolation block for selecting two points of 

10 samples with respect to the signal by the adding means 
and finding the value at the required position from 
linear interpolation, wherein the FIR filter is an FIR 
filter where an impulse response is expressed by a finite 
time length, an impulse response becomes the filter 

15 coefficient, and a transmission function H(z) is 

associated with a transmission function Z(z) of a pre- 
filter, and the filter coefficient is set by performing 
the weighted approximation with respect to a desired 
characteristic in relation to a frequency to be passed 

20 and/or a frequency response of the pre-filter. 

[0031] A ninth aspect of the present invention is an 
audio apparatus including a sampling rate converter, 
wherein the sampling rate converter has a convolution 
processing unit including poly-phase filters able to set 

25 different filter coefficients obtained by poly-phase 
decomposing a predetermined FIR filter and performing 
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convolution processing of input sample signals and a 
poly-phase filter having a selected coefficient, a 
selector for selecting two points of samples required for 
an output sample and selecting the coefficient of the 
5 corresponding poly-phase filter, and a linear 

interpolation block for finding the value at the required 
position from linear interpolation, wherein the FIR 
filter is an FIR filter where an impulse response is 
expressed by a finite time length, the impulse response 

10 becomes the filter coefficient, and a transmission 

function H(z) is related to a transmission function Z (z) 
of a pre-filter, and the filter coefficient is set by 
performing weighted approximation with respect to a 
desired characteristic in relation to a frequency to be 

15 passed and/or a frequency response of the pre-filter. 

[0032] According to the present invention, for example 
a FIR filter is designed. For example, by an initial 
setup, a linear phase FIR filter is set, a band is set, a 
coefficient of the pre-filter is set, any frequency point 

20 to be passed is input, and an initial extreme point is 
set. 

Next, an interpolation polynomial for interpolating 
the amplitude characteristic from a present extreme point 
and a frequency point to be passed is generated. 
25 Next, a new extreme point is determined from the 

amplitude characteristic found from the generated 
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interpolation polynomial . 

These are repeated and it is judged whether or not 
for example the position of the extreme is approximated 
within a desired range. 
5 Then, the filter characteristic is found from the 

approximated amplitude characteristic. 

In this way, in an FIR filter having the 
coefficient set therein, the weighted approximation error 
becomes an equal ripple, and the gain at the pass band is 
10 held at a constant value. 

Further, a designated frequency point can be passed. 
[0033] Then, the up sampler raises the sampling 
frequency U-fold by inserting U-l zero points between 
sample signals . 
15 Next, predetermined convolution processing is 

carried out with respect to the signal multiplied in its 
sampling frequency by U by the convolution processing 
unit including the FIR filter designed as described above. 
Next, it selects two points of samples with respect 
20 to the results of processing of the convolution 

processing unit and finds the value at the required 
position from linear interpolation. 

EFFECT OF THE INVENTION 
[0034] According to the present invention, sampling 
25 rate conversion having a strict cutoff frequency is 
possible. 
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Further, it is possible to avoid chessboard 
distortion . 

Further, any pre-filter can be considered, and it 
is possible to pass any frequency point. 

Further, it is possible to suppress the amount of 
processing to a required lowest limit, and an improvement 
of a processing speed can be achieved. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0035] [FIG. 1] FIG. 1 is a view showing the 
configuration of a transversal type circuit of an FIR 
filter . 

[FIG. 2] FIGS. 2A to 2C are views showing examples 
of a frequency response avoiding chessboard distortion 
and a weighted approximation error in a conventional 
method. 

[FIG. 3] FIGS. 3A to 3C are enlarged views of the 
frequency response and the vicinity of a gain 1 in the 
conventional method . 

[FIG. 4] FIG. 4 is a view showing an example of the 
configuration of a general sampling rate converter. 

[FIG. 5] FIG. 5 is a view of the configuration 
showing a first embodiment of the sampling rate converter 
according to the present invention. 

[FIG. 6] FIG. 6 is a view conceptually showing 
linear interpolation processing of a linear interpolation 
block according to the present embodiment. 

17 
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[FIG. 7] FIG. 7 is a view conceptually showing 
processing for finding a value at a required position by 
the linear interpolation processing of the linear 
interpolation block according to the present embodiment. 

[FIG. 8] FIGS. 8A to 8D are views showing impulse 
responses in four cases where the FIR filter has a linear 
phase . 

[FIG. 9] FIG. 9 is a view showing Q(e 3Q ) and R with 
respect to four cases of the linear phase FIR filter. 

[FIG. 10] FIG. 10 is a view showing an example of 
weighted Chebyshev approximation. 

[FIG. 11] FIG. 11 is a flow chart of a Remez 
Exchange algorithm considering a frequency response of a 
pre-filter according to the present invention. 

[FIG. 12] FIGS. 12A to 12C are views for explaining 
a method for determining a new extreme of a weighted 
approximation error E (e 5 ") . 

[FIG. 13] FIGS. 13A to 13D are views showing a 
frequency response when designating any frequency point 
of the present invention and enlarged views thereof. 

[FIG. 14] FIG. 14 is a view of the configuration 
showing a second embodiment of a sampling rate converter 
according to the present invention. 

[FIG. 15] FIG. 15 is a view of the configuration 
showing a third embodiment of a sampling rate converter 
according to the present invention. 

18 
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[FIG. 16] FIG. 16 is a view of" the configuration 
showing a fourth embodiment of a sampling rate converter 
according to the present invention. 

[FIG. 17] FIG. 17 is a view conceptually showing a 
fifth embodiment of a sampling rate converter according 
to the present invention. 

[FIG. 18] FIG. 18 is a view more specifically 
showing the sampling rate converter conceptually shown in 
FIG. 17. 

[FIG. 19] FIG. 19 is a view for explaining a 
variable used in the linear interpolation block. 

[FIG. 20] FIGS. 10A and 20B are views for 
specifically explaining MasterCount. 

[FIG. 21] FIG. 21 is a view conceptually showing a 
sixth embodiment of a sampling rate converter according 
to the present invention. 

[FIG. 22] FIG. 22 is a view more specifically 
showing the sampling rate converter conceptually shown in 
FIG. 21. 

[FIG. 23] FIG. 23 is a flow chart for explaining a 
first mounting method of a sampling rate converter 
according to the present embodiment. 

[FIG. 24] FIG. 24 is a flow chart for explaining a 
second mounting method of a sampling rate converter 
according to the present embodiment. 

[FIG. 25] FIG. 25 is a block diagram showing an 

19 
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example of* the configuration of an audio apparatus 
employing the sampling rate converter according to the 
present invention. 

DESCRIPTION OF NOTATIONS 
5 [0036] 1... linear phase FIR filter, 2-1 to 2-n- 

1... delay units, 3-1 to 3-n... multipliers, 4... adder, 
h(0) to h(n-l)... filter coefficients, TIN... input 
terminal, TOUT... output terminal, 100, 100A. . . sampling 
rate converter, 101... input terminal, 102... LPF, 103... 

10 up sampler, 104... convolution processing unit, 105... 

linear interpolation block, 106... output terminal, 200, 
200A to 200E... sampling rate converter, 201... input 
terminal, 202... LPF, 203-1 to 203-U. . . convolution 
processing unit, 204-1 to 204-U. . . up sampler, 205-1 to 

15 205-U-l... delay units, 206-1 to 206-U-l... adders, 207... 
down sampler, 208... output terminal, 209... LPF, 210... 
selector, 211... input terminal, 212... FIR filter design 
unit, 213... first memory, 214... LPF convolution 
processing unit, 215... second memory, 216... input 

20 buffer, 217... convolution processing unit, 218... 
counter control unit, 219.... linear interpolation 
processing unit, 220... output terminal, 230... selector, 
231... input terminal, 232... FIR filter design unit, 
233. . . first memory, 234. . . input buffer, 235. . . 

25 convolution processing unit, 236... counter control unit, 
237... linear interpolation processing unit, 238... LPF 
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convolution processing unit, 239. . . second memory/ 24 0. . 
output terminal, 300... audio apparatus, 301... input 
terminal, 302... LR separation circuit (DSB) , 303... 
sampling rate converter (SRC), 304... attenuator (ATT), 
5 305. . . mute circuit (MUTE) , 306. . . output terminal. 

BEST MODE FOR WORKING THE INVENTION 
[0037] Below, a detailed explanation will be given of 
preferred embodiments of the present invention with 
reference to the attached drawings. 
10 [0038] <First embodiment> 

FIG. 5 is a view of the configuration showing a 
first embodiment of a sampling rate converter according 
to the present invention. 

In FIG. 5, "U" and "D" are positive integers 
15 primary with respect to each other, and H(z) indicates a 
transmission function of a FIR filter. Further, the 
upward arrow indicates an up sampler for inserting (U-l) 
zero points between signals. 

Further, Fsi indicates a sampling frequency of 
20 input, and Fso indicates a sampling frequency of output. 
The first embodiment shows an example of the 
configuration where Fso<Fsi . 

[0039] Namely, the sampling rate converter 100 has, as 
shown in FIG* 5, an input terminal 101, low pass filter 
25 (LPF) 102, up sampler 103, convolution processing unit 

104, linear interpolation block 105, and output terminal 
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106. 

[0040] The input terminal 101 receives as input a 
sample signal x(n) having a sampling frequency Fsi. 

[0041] The LPF 102 prevents (suppresses) occurrence of" 
5 folding from the sample signal x(n) of the sampling 

frequency Fsi input from the input terminal 101 since an 
aliasing component is generated and the folding occurs 
where the sampling frequency Fsi of the input is higher 
than the sampling frequency Fso of the output and outputs 
10 the same to the up sampler 103. 

[0042] The up sampler 103 receives the sample data 
x(n) having the sampling frequency Fsi input from the 
input terminal 101 via the LPF 102, inserts 

(interpolates) U-l zero points, raises the sampling 
15 frequency Fsi Up-fold (over sampling) , and outputs the 
sample signal having a sampling frequency UFsi to the 
convolution processing unit 104. 

[0043] The convolution processing unit 104 includes an 
FIR filter designed based on the Remez Exchange algorithm 

20 explained later, performs convolution processing 

(performs band restriction) shown in the following 
equation, and outputs the results of processing to the 
linear interpolation block 105 in the next stage. The 
convolution processing unit 104 interpolates the value by 

25 the convolution of the low pass filter (transmission 
frequency (H(z)) having a cutoff frequency of 1/Up. 

22 
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[0044] [Equation 4] 

y(n') = f; h(i)x(n'-l) - (4) 

1=0 

[0045] Here, h(n) is an impulse response of the FIR 
filter. The output of the convolution (before down 
5 sampler) is the sample obtained by inserting the zero 
point by the up sampler. 

[004 6] The linear interpolation block 105 selects two 
points of samples from the output signal of the 
convolution processing unit 104 having the sampling 
10 frequency UFsi as shown in FIG. 6 and FIG. 7 , finds the 
value at the required position from linear interpolation 
as in the following equation (FIG. 7) , and outputs the 
same as a sampie signal y(m) from the output terminal 106 
[0047] (Equation 5) 
15 y (m) = a, xln A+ jS xln B 

0 =1- a (5) 
[0048] Below, a detailed description will be given of 
the method of design of the FIR filter of the convolution 
processing unit 104 . 
20 [0049] The linear phase FIR filter according to the 

present embodiment can equivalently employ a transversal 
type circuit configuration for example as shown in FIG. 1, 
Note that the filter coefficient h (n) is found from 
the approximated amplitude characteristic by, as will be 
25 explained in detail below, expanding the Remez Exchange 
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algorithm, designating the frequency point to be passed, 
considering the frequency response of the pre-filter, and 
Chebyshev approximating the desired amplitude 
character! s ti c . 
5 [0050] Below, an explanation will be given of a 

specific method of setting the coefficient of the linear 
phase FIR filter according to the present invention in 
sequence with reference to the drawings. 
[0051] As shown in equation (6) , the transmission 
10 function H(z) of the linear phase FIR filter of N taps is 
a filter that stands from a product of the transmission 
function Z(z) of the pre-filter and the transmission 
function K(z) of the equalizer. 
[0052] (Equation 6) 
15 H(z)=Z (z) K(z) (6) 

[0053] Here, it is assumed that the pre-filter and the 
equalizer are linear phase FIR filters of U taps and N- 
(U-l) taps, and the transmission function of the pre- 
filter is previously given. 
20 Further, any Np number of frequency points are 

passed in the frequency domain. Accordingly, the filter 
design of the transmission function H(z) here is to 
determine the transmission function K(z) of the equalizer 
of N-(U-l) taps so that the designated frequency points 
25 are passed and the amplitude characteristic H(e 3G> ) 

approaches the desired amplitude characteristic D (e j£D ) . 
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[0054] The number of taps allocated to the equalizer 
having the transmission function K(z) is set as 3>=N-(U-1). 

The transmission function K(z) of the linear phase 
FIR filter has the linear phases as shown in FIGS. 8A to 
5 8D, so is classified to four cases. 

Specifically/ it is classified into a case 1 of odd 
number taps and even symmetry shown in FIG. 8A, a case 2 
of even number taps and even symmetry shown in FIG. 8B, a 
case 3 of odd number taps and odd symmetry shown in FIG. 
10 8C, and a case 4 of even number taps and odd symmetry 
shown in FIG. 8D . 

[0055] Then, the amplitude characteristic function 
K(e jw ) is held as it is in the case 1, and cases 2 to 4 
are rewritten as follows: 
15 [0056] [Equation 7] 

Casel: ^ a(n)cos(na>) (7-1) 

n-0 



L/2 



Case 2: £ b(n)cos{(n -l/2)cu} 



n-l 

I,/ 2-1 



(7-2) 



cos(— ) ^b(n)cos(nco) 

2 n =0 



(L-l)/2 

Case 3: ^ c(n)sin(nty) 

20 " =1 (7-3) 

(Z.-3)/2 



sin(<y) ]T c(n)cos(r?a>) 



n=0 
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ui ] 
Case 4: ^ d(n)sin{(n -— )co} 

n= ' 2 < 7 " 4 > 

= sin( — )^ d(n)cos(no>) 

2 n =0 



[0057] Namely, the amplitude characteristic 
coefficient K(e j<0 ) is expressed by the product of the 
5 function Q(e 3 °) of a fixed parameter and a cosine series 
P(e 3w ) including the design parameter shown in FIG. 9. 
Below, the upper limit of a sum of Equation (7-1) to 
Equation (7-4) will be expressed by R-l+2xNp. Namely, R 
is computed as in FIG, 9. Further, a(n) ; ~~b(n) ; ~~ c(n) ; 
10 and ~~ d(n) are totally referred to all together as n p(n) f ' 
[0058] When the desired amplitude characteristic is 
defined as D (e jG> ) and the weight with respect to each 
frequency is W(e jw ) , the weighted approximation error is 
defined as follows: 
15 [0059] [Equation 8] 

E (e jtt )=W(e jG> ) {D(e ju )-H(e 3 ") } (8) 
[0060] [Equation 9] 

H(e jw )=K(e 3£,) ) Z (e ja ) 

=Q(e ja> ) P(e j ") Z(e jC) ) (9) 
20 [0061] When entering Equation (9) into Equation (8) , 
the result becomes as follows. 
[0062] [Equation 10] 

E(e j " ) = W(e ;w ){D(e j " ) - P(e j " )} (10) 

[0063] Note that W (e 5w ) and D (e j ") are as follows. 
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[0064] . [Equation 11] 

W (e jo )=W(e jo ) Q(e j<a ) Z(e jw ) (11) 
[0065] [Equation 12] 

D (e j<D )=D(e ju )/Q(e 5 °) Z (e jo ) (12) 
5 [0066] Equation (10) expresses the weighted 

approximation errors of the linear phase FIR filter in 
the four cases of Case 1 to Case 4 . 

The weighted Chebyshev approximation problem is to 
determine such a (n) ; ~~b(n) ; ~~c(n) ; and ~~d(n) of Equation 
10 (7-1) to Equation (7-4) that minimize the maximum value 
of |E(e jG> ) | in the designated frequency band in Equation 
(8) . 

[0067] Below, this will be explained with reference to 
specific examples. 
15 Here, as shown in the following description and FIG. 

10, the amplitude characteristic D (e jw ) is defined. 
[0068] [Equation 13] 

D(e j03 )=l (within error±5i, 0<co<co p ) 
D(e jw )=0 (within error± 5 2 , co s <(d<tz) (13) 
20 [0069] Note that when R is given, the values of 5 1 and 
5 2 cannot be freely designated, but the ratio thereof 
can be designated. 

W(e ju ) is set at a constant value Wl at the pass 
band, while is set at W2 in a blocked frequency and 
25 selected so that Wl 5 1=W2 5 2 stands. For example, it is 
selected so that Wl=l and W2=5l/5 2. At this time, the 
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following reflection theorem stands. 
[0070] <Theorem> 

The required sufficient conditions for the (R-l)th 
order cosine series P(e jG> ) to be the best weighted 
5 Chebyshev approximation with respect to the target 

characteristic in a section (0,7c) of w are: 

(1) E (e 3 ") becomes the extreme at least (R+l) times 

in the section (0,7c). The frequency at which the extreme 

is taken at that time is set as w0<wl<w2<* *<wR-KwR. 
10 (2) The signs of the adjoining extremes are 

different, and the absolute values of all extremes are 

equal. Namely, the following conditions are satisfied. 

[0071] [Equation 14] 

E (e 3<ai ) E (e jwiii )<0 (i=0,l, , R-l) 
15 |E(e jwi ) | = |E(e 3ttiii ) | (i=0,l, ~, R-l) (14) 

[0072] Accordingly, |E(e 3asi )| is equal to the maximum 

value of |E(e jM )| in the section. 

[0073] The technique for obtaining the best Chebyshev 
approximation includes the Remez Exchange algorithm based 

20 on the reflection theorem (refer to Rabiner, L.R. , 

McClellan, J.H. and Parks, T.W.: "FIR Digital Filter 
Design Techniques Using Weighted Chebyshev Approximation'' , 
Proc. IEEE, Vol. 63, April, pp. 595-610, 1975). 

The Remez Exchange algorithm Chebyshev approximates 

25 the desired amplitude characteristic in the frequency 

domain and finds the coefficient of the linear phase FIR 
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filter from the approximated amplitude characteristic. 
[0074] FIG. 11 is a flow chart of the Remez Exchange 
algorithm passing any frequency point and considering the 
frequency response of the pre-filter according to the 
5 present invention. 

A Remez Exchange algorithm considering the 
frequency response of a specific pre-filter becomes as 
follows . 

[0075] <Step 0> 

10 As shown in FIG. 11, first, the initial setup is 

carried out (F101) . In this initial setup, the linear 
phase FIR filter is set, the band is set, the coefficient 
of the pre-filter is set, any frequency point to be 
passed is input, and the initial extreme point is set. 

15 Specifically, the set items are as follows. 

Number of taps, 

Linear phase FIR filter being even symmetric or 
odd symmetric, 

Number of bands, 
20 Frequencies of both ends of each band, 

Desired amplitude value of each band, 
Weighting with respect to each band, 
Coefficient of pre-filter, 

Frequency of point to be passed and amplitude 
25 value (W R+1/ D (e jcoR+1 ) , i=l , Np) , 

Frequency w (0) =w k (0) (k=0, *, R) which becomes the 
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extreme in approximation band. 

Note that the superscript (i) expresses the number. 

of times of repetition. 

[0076] <Step 1> 
5 Next, a Lagrange interpolation polynomial for 

interpolating the amplitude characteristic from the 

present extreme point is generated (F102) . 

The required sufficient condition for the target 

function of the Chebyshev approximation shown in the 
10 above Equation (8) to become the minimum is indicated by 

the reflection theorem. Therefore, based on the 

reflection theorem, the parameter p(n) of the following 

equation is found so that the weighted approximation 

error 5 (i) from the desired amplitude characteristic is 
15 equal at each frequency point and the signs are reflected. 

[0077] [Equation 15] 

p(e j ") = ]T p(n)cos(n<y) (15) 

n=0 

[0078] Namely, the weighted approximation error of 
Equation (9) at the frequency point w U) =Wjc <:l) (K=0 , " , R) 
20 satisfies the following equation: 
[0079] [Equation 16] @ 

»(e jcok(i) ) {D (e 3uk(i) )-P(e j£0k(i) ) }=(-l) k 5 (i) (k=0,l, -, R) 

(16) 

[0080] Below, for simplification, the superscript (i) 
25 is omitted. When modifying Equation (16) , it becomes as 
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follows . 

[Equation 17] 

P(e 3olt ) + (-l) 5c 5 /W<e 3u *)=D <e 3 °*) 

E P(n)cosK) + |^- = D(e^) (k = 0,1,...,R) 

[0082] Then, the equality of the point to be passed in 
the frequency domain is added as a restriction to 
Equation (17) : 

[Equation 18] @ 

P(e joak )=D (e 3c>lc ) 

£ p(n)cos(nty k ) = D(e jaic ), (k = R + l,...,R + Np) (18) 

[0084] When expressing Equation (17) and Equation (18) 
in a matrix, they become as follows. 
[0085] [Equation 19] 
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1 cos(coo) cos(2coo) 

1 cos(cui) cos(2coi) 

• • • 

• • • 

• • • 

1 cos(oji-i) cos(2q>«-i) 

1 cos(COi) cos(2cui) 

1 cos(cutu) COS(2(U*«l) 



cos((R-l)coo) 
cos((R-l)o>i) 



W(e ,w ') 
-1 

W(e Jtl ") 



cos((R-l)co«-i) (-I)"' 1 
W(e Ja " ') 



cos((R-1)o>r) 
cos((R-l)co«4i) 



(-1) R 
W(e'°") 
0 



1 cos(a)mp) cos(2(t)i«p) ••• cos((R-l)coi«p) 0 



P(0) 
P(l) 



p(R-l) 
<5 



D(e'°") 

^ (19) 

D(e JW ') 
D(e ,w ' H ) 

L D(e'"'"') J 

[0086] However, there is a very large amount of 
processing for solving this equation, so 8 will be 
analytically found first. 
[0087] [Equation 20] 

R 



2>,Z>(e)> a ») 



S = 



_ j=0 



[0088] 



X(-iya,/#( e ^) 

7=0 

[Equation 21] 

R 



(20) 



* 1 

a k =fl 1 

[0089] [Equation 22] 

Xj=COS ( CO -j ) 



(21) 



(22) 



32 



S05P0088 

[0090] ak is a remainder factor of elements of the ki- 
th row and (R+l)th column in a matrix F. Note, for W(e 3w ) 
and D (e jw ) , Equation (11) and Equation (12) are used. 

Next, the following equations are set by using this 

5 5 . 

[0091] [Equation 23] 

C k = D(e jQ *)-(-l) k , * 

W(e Jfflk) (23) 

(k = 0,...,R) 
[0092] [Equation 24] 

C k =D (e jwlc ) , (k=R+l, ~, R+Np) (24) 
10 [0093] In order to find the amplitude characteristic 
of the frequency other than the extreme point, the 
Lagrange interpolation polynomial is used this time as 
the interpolation polynomial for interpolation by using 
the extreme point and the frequency point to be passed. 
15 Namely, P(e 3 ") is computed by performing interpolation so 
that the value Ck is taken by wk(k=0, *, R+Np) by using 
the Lagrange interpolation polynomial. 
[0094] [Equation 25] 

R+Np ft 

ZQ( Pk 



20 



p(ey a )= V* } (25) 
Z (- pk 



*=0 * A k 



[0095] [Equation 26] 

A ■ 0 <^> 

j*k 
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[0096] [Equation 27] 

x=cos ( co ) (27) 
[0097] This result corresponds to Equation (19) being 
solved. 
5 [0098] <Step 2> 

A new extreme point is found from the amplitude 
characteristic found from the interpolation polynomial 
(F103) . It is repeatedly decided whether or not the 
optimum approximation is obtained (F104) . 
10 Each extreme point wk of the result of above step 1 

does not always become the extreme of the weighted error 
function E (e jw ) . A point where |E(e jw ) I becomes larger 
than 5 (1> sometimes exists. Therefore, the new extreme 
point w u+1) is determined by the all point simultaneous 
15 replacement method. 

<A11 point simultaneous replacement method> 
Based on the following equation, the extreme of the 
weighted approximation error computed from the extreme 
point used for interpolation is searched for over the 
20 approximation band. This is defined as a new extreme 

point w (i+1) =w k (i+1) (k=0,l, R) , then the routine returns 

to the processing of step 1 . 
[0099] [Equation 28] 

E (e ju ) =W (e j °) { D (e 3w ) -P (e 3 ") } (28) 
25 [0100] It is assumed that the optimum approximation is 
obtained when the position of the extreme no longer 
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changes. This is the end condition of the repetition. The 
routine then proceeds to the processing of the next step 
3. 

[0101] FIGS. 12A to 12C are conceptual views of the 
5 all point replacement method. 

Simply explaining this, the black dots in FIGS. 12A 
to 12C express extreme points used for the interpolation. 
The weighted approximation error E (e jw ) found from this 
extreme point corresponds to the solid line. 

10 As shown in FIG. 12A, the value of the weighted 

approximation error at the extreme point of a black dot 
becomes a white dot, but the actual extreme is the 
frequency indicated by a square. Therefore, the routine 
returns to the processing of step 1 while defining the 

15 frequency indicated by the square as the new extreme 
point. 

Furtheir, as shown in FIG. 12B, the frequency of the 
extreme point used for the interpolation and the 
frequency of the actual extreme are offset, therefore the 
20 frequency indicated by the square is defined as the new 

extreme point, then the routine returns to the processing 
of step 1 . 

Then, as shown in FIG. 12C, the repetition ends 
when the extreme point used for the interpolation and the 
25 extreme point (white dot) of the actual weighted 
approximation error become the same. 
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[0102] <Step 3> 

The coefficient of the linear phase FIR filter is 
found from the approximated amplitude characteristic 
(F105) . 

5 When finding the impulse response h(n) of N taps 

from the optimum approximate function P(e jG) ), it is found 
from the following equation in place of finding the same 
from p (n) . 

[0103] [Equation 29] 
10 H(e jtD )=P(e 3 ") Q(e 3u ) Z (e jtt ) (29) 

[0104] [Equation 30] 
Case 1 : 

N-1 

h(n) = ^{H(0)} + 2£ (-l) k H(^k)cos(— k(n + -))} (30) 

[0105] [Equation 31] 
15 Case 2 : 

9 2 ?7T 1 ?7T 1 1 

h(n) = — V (-l) k H(— (k + -))sin(— (k + -)(n + -)) (31) 

[0106] [Equation 32] 
Case 3 : 

N-1 

O 2 977- 777- 1 

h( n ) = -— V (-l) k H(— k)sin(— k(n + -)) (32) 

20 [0107] [Equation 33] 
Case 4 : 



h(n)= |i (-l) k H(^(k + i))cos(^(k + |)(n + i)) (33) 
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[0108] Further, when finding the impulse response k(n) 
of the equalizer of L=N-(U-1) taps, it is computed from 
the following equation. 
[0109] [Equation 34] 
5 K (e jtt )=P(e j ") Q(e jo ) (34) 

[0110] [Equation 35] 
Case 1 : 

L-l 

k(n) = -{K(0)} + 2X (-l) k K(-^k)cos(— k(n + -))} (35) 

[0111] [Equation 36] 
10 Case 2 : 

9 2 9^ i 9^ i i 

k ( n ) = f Z (-l) k K(^(k + i-))sin(^(k + -)(n + -)) (36) 

[0112] [Equation 37] 
Case 3 : 

kl 

k ( n ) = ~f£ (-l) k K(^k)sin(^k(n + |)) (37) 

15 [0113] [Equation 38] 
Case 4 : 

9 2 9^ i 9^ i i 

k ( n ) = f E (-l) k K(^(k + ^))cos(^(k + ^)(n + i)) (38) 

[0114] If the transmission function Z(z) of the pre- 
filter is 1 as shown in the following equation, it is the 
20 same as the Remez Exchange algorithm passing through any 
frequency point. 
[0115] (Equation 39) 
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Z(z)=l (39) 
[0116] Further, in the case where there is no 
frequency point to be passed through the pre-filter, when 
Np=0, it is the same as a Remez Exchange algorithm 
5 considering the frequency response of the pre-filter. 

[0117] Further, when the transmission function Z(z) of 
the pre-filter is 1 as shown in the following equation 
and there is no frequency point to be passed, when Np=0, 
it is the same as an ordinary Remez Exchange algorithm. 
10 [0118] (Equation 40) 

Z(z)=l (40) 
[0119] FIGS. 13A to 13D are views showing the 
frequency response of the low pass filter designed by the 
Remez Exchange algorithm expanded so that any frequency 
15 point is passed and the frequency response of the pre- 
filter can be considered for the following specifications. 

Note that in the following explanation, the zero 
point for avoiding chessboard distortion will be handled 
as a pre-filter. 
20 The frequency response of the pre-filter is 

expressed as follows. 
[0120] [Equation 41] 

Z(e 3 ") = (l+e- jft, +e- 23a> + " +e" j< " ^ ) 



(C/-l)/2 

1+ 2]2cos(/77cy), U\odd 

U/2 

2 cos{(w + 1 / 2)o))}, U : even 
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[0121] The specifications will be shown below. 
[0122] <Linear phase FIR filter> 
24 taps 
5 Even symmetric 

Pre-filter U=3 of Equation (2) (adjusted so that 
DC gain becomes U) 

<Design method> 
Designed by Remez Exchange algorithm passing 
10 through any frequency point and considering frequency 
response of pre-filter. 
[0123] [Table 1] 
Band 



Band 


Frequency range 


Gain 


Weight 


1 


0 = CO = 0 . 3 7E 


3 


1 


2 


0 . S TZ = CO = 71 


0 


1 



15 [0124] [Table 2] 

Designated frequency point 



Frequency 


Gain 


co =3 


3 



[0125] FIG. 13A is a view showing the frequency 

20 response displayed by decibels, FIG. 13B is a view 

showing the frequency response displayed while keeping 
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the value as it is, FIG. 13C is a view enlarging the 
vicinity of the gain 3, and FIG. 13D is a view enlarging 
the vicinity of the gain 0. 

In FIGS. 13A to 13D, the dotted lines indicate the 
5 frequency characteristics of the pre-filter and the 
frequency characteristics of the equalizer, the solid 
curves indicate the finally obtained frequency 
characteristics (Proposed H(z)), the vertical solid lines 
indicate the frequencies (zero point) where H(z)=0 must 
10 stand so as to avoid chessboard distortion, and the black 
dots indicate sections of the band. 

[0126] From FIG. 13A, it can be confirmed that the 
gain of the pass band maintains a constant value and 
passes the zero point for avoiding chessboard distortion. 
15 Further, it can be confirmed from FIG. 13C that the 

designated frequency point is passed. 

Further, it can be confirmed from FIGS. 13C and 13D 
that an equal ripple is held. 

[0127] Namely, a low pass filter designed by a Remez 
20 Exchange algorithm which passes any frequency point and 
is enlarged so that the frequency response of the pre- 
filter can be considered can give a good frequency 
response characteristic . 

[0128] Next, an explanation will be given of the 
25 operation of the sampling rate converter having the above 
configuration . 
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[0129] The sample data x(n) having the sampling 
frequency Fs input from the input terminal 101 is input 
to the LPF 102. 

In the LPF 102, the occurrence of folding from the 
5 sample signal x(n) of the sampling frequency Fsi input 

from the input terminal 101 is prevented (suppressed) and 
the signal is output to the up sampler 103. 

In the up sampler 103, (U-l) zero points are 
inserted between signals, the sampling frequency Fsi is 
10 raised U-fold, and the sample signal having the sampling 
frequency UFsi is output to the convolution processing 
unit 104. 

The convolution processing unit 104 performs the 
convolution processing based on Equation (3) , restricts 

15 the band of the sample signal, and supplies the same to 
the linear interpolation block 105 in the next stage. 

At this time, the output of the convolution (before 
down sampler) is the sample obtained by inserting the 
zero point by the up sampler. 

20 Then, the linear interpolation block 105 selects 

two points of samples from the output signal of the 
convolution processing unit 104 of the sampling frequency 
UFsi. Then, the value at the required position is found 
from linear interpolation. 

25 By this, the sample signal y (m) having the 

sampling frequency Fso is output from the output terminal 
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106. 

[0130] In the sampling rate converter 100 having the 
above function, the FIR filter of the convolution 
processing unit 104 is a linear phase FIR filter in which 
5 the impulse response becomes the filter coefficient, the 
transmission function H(z) is associated with the 
transmission function Z (z) of the pre-f ilter and the 
transmission function K(z) of the equalizer and 
associated with the frequency point to be passed and the 

10 frequency response of the pre-f ilter, and the filter 

coefficient is set based on the amplitude characteristic 
of the equalizer obtained by performing weighted 
approximation with respect to the desired characteristic 
by using a Remez Exchange algorithm passing through any 

15 frequency point and considering the frequency response of 
the pre-filter with reference to the frequency point to 
be passed and the frequency response of the pre-filter, 
therefore the present sampling rate converter has the 
following advantages. 

20 Namely, chessboard distortion can be avoided. 

Further, any pre-filter can be considered, and any 
frequency point can be passed. 

[0131] Further, provision was made of the up sampler 
103 for inserting (U-l) zero points between signals and 
25 raising a sampling frequency Fsi U-fold, the convolution 
processing unit 104 for interpolating the value by the 
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convolution of the low pass filter (transmission function 
H(z)) having the cutoff frequency of 1/Up, and the linear 
interpolation block 105 for selecting two points of 
samples from the output signal of the convolution 
5 processing unit 104 having the sampling frequency UFsi 
and finding the value at the required position from the 
linear interpolation, therefore, as shown in Table 3, 
there is the advantage that sampling rate conversion 
wherein the cutoff frequency is severe becomes possible. 
10 [0132] Table 3 shows the frequency band supported by 
the sampling rate converter according to the present 
embodiment . 
[0133] [Table 3] 

Table 3. Frequency Supported by Sampling Rate Conversion 
15 of Present Invention 





Output (Fso) 


Input (Fsi) 


8 kHz 


11 .025 
kHz 


12 kHz 


22.05 
kHz 


24 kHz 


32 kHz 


44.1 
kHz 


48 kHz 


8 kHz 




O 


A 


O 


A 


A 


O 


A 


11.025 kHz 


O 




O 


. A 


O 


O 


A 


O 


12 kHz 


A 


O 




O 


A 


A 


O 


A 


22.05 kHz 


O 


A 


O 




O 


o 


A 


O 


24 kHz 


A 


O 


A 


O 




A 


O 


A 


32 kHz 


A 


O 


A 


O 


A 




O 


A 
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44.1 kHz 






o 




o 


o 




o 


48 kHz 


A 


O 


A 


O 


A 


A 


O 




96 kHz 


A 


O 


A 


O 


A 


A 


O 


A 



[0134] In Table 3, the sampling frequencies Fsi of the 
inputs are 8 kHz, 11.025 kHz, 12 kHz, 22.05 kHz, 24 kHz, 
32 kHz, 44.1 kHz, 48 kHz, and 96 kHz, and the sampling 
5 frequencies Fso of the outputs are 8 kHz, 11.025 kHz, 12 
kHz, 22.05 kHz, 24 kHz, 32 kHz, 44.1 kHz, and 48 kHz. 
[0135] In Table 3, the portions of "a" are frequencies 
which can be achieved by the sampling rate converter 
based on the general multi-rate signal processing. In 

10 general, those based on multi-rate signal processing have 
better characteristic of attenuation and are easier in 
control of the pass band etc. , but achievement is 
possible also by the sampling rate converter according to 
the present embodiment. 

15 [0136] Further, according to the first embodiment, the 
LPF 102 is provided between the input terminal 101 and 
the up sampler 103, therefore, when the sampling 
frequency Fsi of the input is higher than the sampling 
frequency Fso of the output, the occurrence of the 

20 aliasing component and the occurrence of folding can be 
suppressed. 

[0137] <Second Embodiment> 

FIG. 14 is a view of the configuration showing a 
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second embodiment of the sampling rate converter 
according to the present invention. 

[0138] The difference of the second embodiment from 
the first embodiment explained above resides in that an 
5 LPF 107 is provided between the linear interpolation 
block 105 and the output terminal 106 in place of 
providing the LPF 102 between the input terminal 101 and 
the up sampler 103, and the configuration is made so that 
an imaging component occurs when the sampling frequency 
10 Fsi of the input is lower than the sampling frequency Fso 
of the output, and occurrence of a frequency component 
which does not exist in the original data is prevented 

(suppressed) . 

[0139] The rest of the configuration is the same as 
15 that of the first embodiment explained above. 

[0140] According to the second embodiment, there are 
the advantages that the occurrence of a frequency 
component which does not exist in the original data can 
be prevented (suppressed) , and, in the same way as the 
20 first embodiment, chessboard distortion can be avoided, 
and the optional pre-filter can be considered, any 
frequency point can be passed, and sampling rate 
conversion wherein the cutoff frequency is severe becomes 
possible . 

25 [0141] Note that it is also possible to provide the 

LPF 102 between the input terminal 101 and the up sampler 
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103 and provide the LPF 107 between the linear 
interpolation block 105 and the output terminal 106. 
[0142] <Third Embodiment 

FIG. 15 is a view of the configuration showing a 
5 third embodiment of the sampling rate converter according 
to the present invention. 

[0143] Note that, in FIG. 15, "U" and "D" are positive 
integers primary with respect to each other, and R(z) 
indicates the transmission function of the poly-phase 
10 filter. Further, the upward arrow indicates the up 

sampler inserting (U-l) zero points between signals. 
[0144] The difference of the third embodiment from the 
first embodiment explained above resides in that the FIR 
filter is so-called "poly-phase decomposed" to U number 
15 of filters and a poly-phase configuration is employed. 

[0145] Here, an explanation will be given of the poly- 
phase decomposition first. 
[014 6] <Poly-phase decomposi tion> 

The sampling rate converter restricts the band by 
20 using the (N-l)th order FIR filter (transmission function 
H(z)) and interpolates the sample of the zero point 
portion. 

[0147] [Equation 42] 

#( z ) = £/,(„) z -« (42) 

25 [0148] The cutoff frequency co c of this FIR filter 
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becomes as follows: 
[0149] [Equation 43] 

wc= 7i /maxU , D (43) 
[0150] The configuration of FIG. 5 can be equivalently 
5 expressed as in FIG. 15 by the poly-phase configuration. 
The transmission function H(z) of the filter of Equation 
(42) and the transmission function Ri (z) of the poly- 
phase filter are expressed by the following relationship. 
[0151] [Equation 44] 
10 H(z)=R u - 1 (z u )+z- 1 R u - 2 (z u )+ +z- (U " 1) R 0 (z u ) 

= ^RKz u )z^ (44) 

n=0 

[0152] [Equation 45] 

N/U-l 

Ri(z)= 2 h(Un-i-l + U)z" n (45) 

n=0 

[0153] Note that the calculation is carried out by 
15 raising N to a whole multiple of U or, where insufficient 

for NU taps, assuming that the coefficient of 0 exists. 

[0154] The input/output relationship of the 

interpolator shown in FIG. 15 becomes as follows by the 

convolution and the processing of the up sampler when 
20 expressing the impulse response of the poly-phase filter 

Ri (z) as ri (n) . 

[0155] [Equation 46] 
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y(m) = 



[0156] Note that k is an integer, and xi (m) becomes as 
follows : 
5 [0157] [Equation 47] 

co 

xi(m) = ^ n(n)x(m - h) (47) 

[0158] The sampling rate converter 200 according to 
the third embodiment has, as shown in FIG. 15, an input 
terminal 201, LPF 202, convolution processing units 203-1 

10 to 203-U, up samplers 204-1 to 204-U, delay units 205-1 
to 205-U-l, adders 206-1 to 206-U-l, a linear 
interpolation block 207, and output terminal 208. 

Among these components, the adding means is 
configured by the delay units 205-1 to 205-U-l and the 

15 adders 206-1 to 206-U-l. 

[0159] The input terminal 201 receives as input the 
sample signal x(n) having the sampling frequency Fs . 
[0160] When the sampling frequency Fsi of the input is 
higher than the sampling frequency Fso of the output, an 

20 aliasing component is generated and folding occurs, 

therefore the LPF 202 prevents (suppresses) the occurrence 
of folding from the sample signal x(n) of the sampling 



m 

*t/-i(— X 



m = kU + (U-\) 

m = kU + l 
m = kU 



(46) 
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frequency Fsi input from the input terminal 201 and 
outputs the results to the convolution processing units 

203- 1 to 203-U. 

[0161] The convolution processing units 203-1 to 203-U 
5 include the poly-phase filters obtained by poly-phase 
decomposing the FIR filter designed based on the above- 
explained Remez Exchange algorithm, perform the 
convolution processing between the input sample signals 
and poly-phase filters based on the above Equation (47) , 
10 and output the results of processing to the up samplers 

204- 1 to 204 -U in the next stage. 

[0162] The up samplers 204-1 to 204 -U receive the 
output sample signals of the convolution processing units 
203-1 to 203-U, insert (U-l) zero points, raise the 

15 sampling frequencies Fsi U-fold, and output sample 
signals having the sampling frequency UFsi . 
[0163] The delay unit 205-1 delays the sample signal 
having the sampling frequency UFs by the up sampler 204-1 
by a predetermined time and outputs the result to the 

20 adder 206-1. 

Further, the delay units 205-2 to 205-U-l delay the 
output signals of the adders 206-1 to 206-U-2 by the 
predetermined time and output the results to the adders 
206-2 to 206-U-l. 

25 [0164] The adder 206-1 adds the sample signal having 
the sampling frequency UFsi by the up sampler 204-2 and 
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the sample signal delayed by the delay unit 205-1 and 
outputs the result to the delay unit 205-2. 

The adder 206-2 adds the sample signal having the 
sampling frequency UFsi by the up sampler 204-3 and the 
5 sample signal delayed by the delay unit 205-2 and outputs 
the result to the delay unit 205-3. 

In the same way, the adder 206-U-l adds the sample 
signal having the sampling frequency UFsi by the up 
sampler 204 -U and the sample signal delayed by the delay 
10 unit 205-U-l and outputs the result to the linear 
interpolation block 207. 

[0165] The linear interpolation block 207 selects two 
points of samples from the output signal of the adder 
206-U-l, that is, the signal obtained by adding sample 

15 signals having sampling frequencies multiplied by U by 

the up samplers 204-1 to 204-U-l, finds the value at the 
required position (FIG. 7) from the linear interpolation 
as in above Equation (5) , and outputs the result as the 
sample signal y(m) from the output terminal 208. 

20 [0166] In the sampling rate converter 200, sample 

signals x(n> having frequencies Fs input from the input 
terminal 201 are input in parallel to the convolution 
processing units 203-1 to 202-U including the poly-phase 
filters after the occurrence of folding is prevented 

25 (suppressed) at the LPF 202. 

[0167] The convolution processing units 203-1 to 203-U 
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perform the convolution processing of the input sample 
signals and poly-phase filters and supply the results of 
processing to the up samplers 204-1 to 204-U in the next 
stage . 

5 The up samplers 204-1 to 204-U insert (U-l) zero 

points between output sample signals of the convolution 
processing units 203-1 to 203-U, raise the sampling 
frequency Fsi U-fold, and output sample signals having 
the sampling frequency UFsi. 

10 The output signals of the up samplers 204-1 to 204- 

U are delayed by the delay units 205-1 to 205-U-l and 
adders 206-1 to 206-U-l, then cumulatively added and 
supplied to the linear interpolation block 207. 

Then, the linear interpolation block 207 selects 

15 two points of samples from the output signal of the adder 
206-U-l having the sampling frequency UFsi. Then, the 
value at the required position is found from the linear 
interpolation. 

By this, the sample signal y(m) having the sampling 

20 frequency Fso is output from the output terminal 208. 

[0168] According to the third embodiment, in addition 
to the effects of the first embodiment explained above, 
there are the advantages that the amount of processing 
can be suppressed to the required lowest limit, and an 

25 improvement of the processing speed can be achieved. 
[0169] <Fourth Embodiment 
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FIG. 16 is a view of the configuration showing a 
fourth embodiment of the sampling rate converter 
according to the present invention. 

[0170] The difference of the fourth embodiment from 
5 the above-explained third embodiment resides in that an 
LPF 209 is provided between the linear interpolation 
block 207 and the output terminal 208 in place of 
providing the LPF 202 between the input terminal 201 and 
the convolution processing units 203-1 to 203-U, and the 

10 configuration is made so that generation of an imaging 
component where the sampling frequency Fsi of the input 
is lower than the sampling frequency Fso of the output 
and generation of a frequency component which does not 
exist in the original data are prevented (suppressed) . 

15 [0171] The rest of the configuration is the same as 
that of the third embodiment. 

[0172] According to the fourth embodiment, in addition 
to the effects of the second embodiment explained above, 
there are the advantages that the amount of processing 

20 can be suppressed to the required lowest limit and an 
improvement of the processing speed can be achieved. 
[0173] Note that it is also possible to provide the 
LPF 202 between the input terminal 201 and the 
convolution processing units 203-1 to 203-U and provide 

25 the LPF 209 between the linear interpolation block 207 
and the output terminal 208. 
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[0174] <Fifth Embodiments 

FIG. 17 is a view conceptually showing a fifth 
embodiment of the sampling rate converter according to 
the present invention. 
5 [0175] The difference of the fifth embodiment from the 
above -explained third embodiment resides in that the 
position of the sample required for the linear 
interpolation can be specified in the linear 
interpolation block, therefore, the configuration is made 

10 so that the linear interpolation block is regarded as a 
selector 210, the poly-phase filter coefficient set 
corresponding to the sample output by the selector 210 is 
selected, and the convolution processing of the poly- 
phases based on the above Equation (4 7) is carried out. 

15 [0176] In this way, it is possible to regard the 

linear interpolation block as the selector 210, select 
the poly-phase filter corresponding to the output sample, 
and carry out only the convolution processing of the 
poly-phases of Equation (47) . 

20 By doing this, it becomes unnecessary to perform 

unnecessary computation. 

[0177] The configuration of the sampling rate 
converter 200B in this case basically has, as shown in 
FIG. 17, an input terminal 201 to which a sample signal 
25 x(n) having a sampling frequency Fs is input, an LPF 202 
for suppressing the generation of the aliasing component 
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and the occurrence of folding when the sampling frequency 
Fsi of the input is higher than the sampling frequency 
Fso of the output, convolution processing units 203 (-1 
to -n) for performing the convolution processing 
5 (Equation (47)) of the input samples and the poly-phase 

filters obtained by decomposing the filter to poly-phases, 
a selector 210 for selecting the poly-phase filter 
corresponding to the output sample, and an output 
terminal 208 for outputting the sample signal y(m) having 
10 the frequency converted to the desired sampling frequency. 
[0178] FIG. 18 is a view more specifically showing the 
sampling rate converter 200B conceptually shown in FIG. 
17. 

[0179] This sampling rate converter 200C has, as shown 
15 in FIG. 18, an input terminal 211 to which the sample 

signal x(n) having the sampling frequency Fs is input, an 
FIR filter design unit 212 having a coefficient set based 
on the algorithm of FIG. 11 explained above, a first 
memory 213 for storing the coefficient set of the poly- 
20 phase filters by the FIR filter design unit 212, a 
convolution processing unit 214 for performing the 
convolution processing of the LPF for suppressing the 
generation of the aliasing component and the occurrence 
of folding when the sampling frequency Fsi of the input 
25 is higher than the sampling frequency Fso of the output, 
a second memory 215 for storing the LPF coefficient, an 
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input buffer 216 for storing input data of the 
convolution processing unit 214, a convolution processing 
unit 217 for performing the convolution processing of the 
input data and the poly-phase filter decomposed to the 
5 poly-phases based on the input data stored in the input 
buffer 216 and the coefficient read out from the first 
coefficient memory 213, a counter control unit 218 having 
a role for selecting two points of samples required for 
the output sample and selecting the corresponding poly- 

10 phase filter and finding the coefficient of the linear 

interpolation, a linear interpolation processing unit 219 
for finding the value at the required position from the 
output of the convolution processing unit 217 based on 
the coefficient of the linear interpolation by the 

15 counter control unit 218, and an output terminal 220 for 
outputting the sample signal converted in frequency to 
the desired sampling frequency. 

[0180] The counter control unit 218 is included in the 
linear interpolation block and has a variable master 
20 counter (MasterCounter) . 

[0181] FIG. 19 is a view for explaining the variables 
used in the linear interpolation block. 

[0182] As shown in FIG. 19, the linear interpolation 
block includes, as variables, a MasterCount as the master 
25 counter for determining the operation of the linear 

interpolation, a CountlOOO as the counter of the decimal 
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point or less to be added to the master counter, an added 
value of an integer portion to be added to the master 
counter, and a CountAmari as the added value of errors of 
the decimal point or less to be added to CountlOOO. 
5 [0183] FIGS. 20A and 20B are views for specifically 
explaining the MasterCount. 

[0184] MasterCount is comprised of at least 26 bits as 
shown in FIGS. 20A and 20B, in which the 16 bits from the 
0-th bit to 15-th bit indicate the coefficient of the 
10 linear interpolation, the 7 bits from the 16-th bit to 
the 22nd bit indicate the number of the poly-phase 
coefficient set, and the 3 bits from the 23rd bit to the 
25-th bit indicate the sample number to be read and 
fetched. 

15 [0185] <Sixth Embodiment 

FIG. 21 is a view conceptually showing a sixth 
embodiment of the sampling rate converter according to 
the present invention. 

[0186] The difference of the sixth embodiment from the 
20 above fourth embodiment resides in the configuration that 
the position of the sample required for the linear 
interpolation can be specified in the linear 
interpolation block, therefore the linear interpolation 
block is regarded as the selector 230, the poly-phase 
25 filter coefficient set corresponding to the sample output 
by the selector 230 is selected, and the convolution 
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processing of the poly-phases based on the above Equation 
(47) is carried out. 

[0187] In this way, only the convolution processing of 
Equation (4 7) may be carried out by regarding the linear 
5 interpolation block as the selector 230 and selecting the 
poly-phase filter corresponding to the output sample. 

By doing this, it becomes unnecessary to perform 
unnecessary computation . 

[0188] The configuration of the sampling rate 

10 converter 200D in this case basically has, as shown in 
FIG. 21, an input terminal 201 to which a sample signal 
x(n) having a sampling frequency Fs is input, convolution 
processing units 203 (-1 to -n) for performing the 
convolution processing (Equation (47) ) of the input 

15 samples and the poly-phase filter decomposed to poly- 
phases, a selector 230 for selecting the poly-phase 
filter corresponding to the output sample, an LPF 209 for 
suppressing generation of the imaging component and 
occurrence of a frequency component not existing in the 

20 original data when the sampling frequency Fsi of the 
input is lower than the sampling frequency Fso of the 
output, and an output terminal 208 for outputting a 
sample signal y (m) having a frequency converted to the 
desired sampling frequency. 

25 [0189] FIG. 22 is a view more specifically showing the 
sampling rate converter 200D conceptually shown in FIG. 
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21. 

[0190] This sampling rate converter 200E has, as shown 
in FIG. 18, an input terminal 231 to which a sample 
signal x(n) having a sampling frequency Fs is input, an 
5 FIR filter design unit 232 having a coefficient set based 
on the algorithm of FIG. 11 explained above, a first 
memory 233 for storing the coefficient set of the poly- 
phase filters by the FIR filter design unit 232, an input 
buffer 234 for storing the input data of the convolution 

10 processing unit 235, a convolution processing unit 235 
for performing the convolution processing of the input 
data and the poly-phase filter decomposed to poly-phases 
based on the input data stored in the input buffer 234 
and the coefficient read out from the first coefficient 

15 memory 233, a counter control unit 236 having a role of 
selecting two points of samples required for the output 
sample and selecting the corresponding poly-phase filter 
and finding the coefficient of the linear interpolation, 
a linear interpolation processing unit 237 for finding 

20 the value at the required position from the output of the 
convolution processing unit 235 based on the coefficient 
of the linear interpolation by the counter control unit 
236, a convolution processing unit 238 for performing the 
convolution processing of the LPF for suppressing the 

25 generation of an imaging component and the generation of 
a frequency component not existing in the original data 
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when the sampling frequency Fsi of the input is lower 
than the sampling frequency Fso of the output, a second 
memory 239 for storing the LPF coefficient, and an output 
terminal 240 for outputting the sample signal converted 
5 in frequency to the desired sampling frequency. 

[0191] Note that the concrete nature of the counter 
control unit 236 is the same as that explained with 
reference to FIG. 19, and FIGS. 20A and 20B. 
[0192] Below, an explanation will be given of the 

10 mounting method and a specific example of the sampling 

rate converter according to the embodiments in sequence. 

Here, an explanation will be given of the mounting 
method and a specific example of the sampling rate 
converter of FIG. 5, FIG. 17, and FIG. 18 and the 

15 mounting method and a specific example of the sampling 
rate converter of FIG. 14, FIG. 21, and FIG. 22 in 
sequence . 

[0193] <Mounting Method of Sampling Rate Converter 
(FIG. 5, FIG. 17, FIG. 18)> 

20 When realizing the sampling rate converter shown in 

FIG. 5, it is realized by using the poly-phase 
configuration shown in FIG. 17. In the method of 
realization explained below, the required smallest 
computation is satisfactory by selecting the required 

25 input with respect to the output and selecting the poly- 
phase filter . 
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[0194] FIG. 23 is a flow chart for explaining the 
mounting method of the sampling rate converter according 
to the present embodiment. 

The specific processing becomes as follows. Note 
that, for convenience of the explanation, C language-like 
expressions are used, but does not matter if the mounting 
means is hardware, software, or a mixture of hardware and 
software . 

[0195] <Step 10> 

As shown in FIG. 23, first, the initial setup is 
carried out (F201) . In this initial setup, the number of 
input data, the number of output data, the center of the 
input data, the number of taps of the poly-phase filter, 
and the number of the coefficient set of the poly-phase 
filters are set. Specifically, the initial value of the 
counter control unit 215 is set. 

They are specifically as follows. 
Input data number: InputSample; 
* Number of taps : Tap ; 

Center of input data: InputOf f set=0 ; 

Number of taps of poly-phase filter: 
PolyTap= (Tap+Up-1) /Up; 

Number of coefficient set of poly-phase filter: 
Coeff Count; and 

Linear interpolation block: 

The method of finding CountAdd and CountAmari is 
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defined from the sampling frequency Fsi of the input and 
the sampling frequency Fso of the output. 

Only integer portion*/ 

CountAdd=(int) (Fso/Fsi*2~23) ; 
5 1000 ties of decimal point or less*/ 

CountAmari=( (Fso/Fsi*2~23) -CountAdd) *1000; 
[0196] When the frequency is converted from 44.1 kHz 
to 48 kHz, it specifically becomes as follows: 

Example: Convert 44.1 kHz — > 48 kHz*/ 
10 CountAdd=9130457; 

CountAmari=687 ; 

The initial value of CountlOOO is determined as 
follows . 

Countl000=500;/*for rounding*/ 
15 [0197] <Step 11> 

Next, the low pass FIR filter is designed by the 
Remez Exchange algorithm already explained with reference 
to for example FIG. 11 (F202) . 

Here, a detailed description thereof will be 
20 omitted. 

[0198] <Step 12> 

Next, the poly-phase filter is prepared (F203) . 
Namely, from the transmission function H(z) of the FIR 
filter, the transmission function Ri(z) of the poly-phase 
25 filter is found by using a development equation of the 
above Equation (42) . 
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As in the following equation, a normalization 
coefficient RiNormai of each poly-phase is found. 
[0199] [Equation 48] 

PofyTap-l 

E Coe /(*0 (48) 
[0200] Then, coefficient inversion for convolution is 
carried out. Namely, in order to perform the convolution 
of above Equation (4 7) , the coefficients of the poly- 
phase filters are rearranged to an inverse sequence. 
[0201] Coefficient set of first poly-phases 

So that the position of the center of taps becomes 

the first output*/ 

It is necessary to select the set of coefficients. 

*/ 

Coef fCount= (PolyTap*/Up/2) %Up; 

Set initial value of MasterCount. */ 

Number of coefficient set of poly-phases including 
the center value of MasterCount=H (z) 
Coef f Count) «16 : 

Note that, «16 indicates a 16-bit shift. 
[0202] <Step 13> 

Next, the poly-phase filter and the center of input 
are determined (F204) . 

The coefficient LinearCoef of the linear 
interpolation and the number CoefCount of the coefficient 
set of the poly-phases are found as follows: 
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/*Coef f icient set and linear interpolation 
coefficient are found from MasterCount*/ 

CoefCount=( (MasterCount&0x007f f f f f ) »16) &0x0000007f 



LinearCoef=( (MasterCount&0x007f f f f f ) SOxOOOOffff ) ; 
[0203] <Step 14> 

Filter processing by the low pass filter LPF is 
carried out (F209) * 

Specifically, after the convolution is carried out 
as in Equation (49) and rounded as in Equation (50) next, 
it is normalized as in Equation (51) , then the clipping 
is carried out as in Equation (52) . 
[0204] [Equation 49] 



Convolution 



LPFTap-l 



Yconv = £ LPF(k)xx(InputOfFset) + k-LPFT ap /2 



(49) 



k=0 



[0205] 



[Equation 50] 



Rounding off 




Yconv < 0 



Yconv > 0 



(50) 



[0206] 



[Equation 51] 



Normalization 




(51) 



[0207] 



[Equation 52] 



Clipping 
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Src(n) = < 



CLIP_MAX 
CLIP MIN 



Nonnal 



Y Nonnal > CLIMAX 
Y Nonnal <CLIP_MIN 
else 



(52) 



[0208] In the case of 16 bits, CLIP MAX=32767 and CLIP 
MIN=32768 stand. 

In the case of 24 bits, CLIP MAX=8388607 and CLIP 
5 MIN=8388608 stand. 

In the case of signed 8 bits, CLIP MAX=127 and CLIP 
MIN=128 stand. 

In the case of unsigned 8 bits, CLIP MAX=255 and 
CLIP MIN=0 stand. 
10 [0209] <Step 15> 

Next, the convolution of the above Equation (47) is 
carried out (F205) . 

Two points of samples required in the linear 
interpolation block are found. In the computation of the 
15 convolution, numbers of the coefficient set of poly- 
phases and the center of the input data are computed by 
using the following values for a left sample InA and a 
right sample InB as shown in Table 4 . 
[0210] [Table 4] 



Output 


Number of poly -phase 
coefficient set 


Center of input 
data 


Remarks 


InA 


Coef Count 


InputOf f set 




InB 


(Coef Count+l ) 
0 


InputOf f set 
InputOffset+1 


Time of 

( Coef Count+l )<1 2 8 
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Time of 

(Coef Count+1 ) =1 28 



[0211] Specifically, after the convolution is carried 
out as in Equation (53) and the result rounded off as in 
Equation (54) , the result is normalized as in Equation 
5 (55) , then clipped as in Equation (56) . 
[0212] [Equation 53] 
Convolution 

PolyTap-l 

In conv = 2 CoefCount(k)xSrc(InputOflfset + (k - PolyTap/2)) (53) 

k=0 

[0213] [Equation 54] 

10 Rounding off 

In conv +(R Normal /2) 

conv — 



Konv-CRNo^a,^) < 0 (54) 



[0214] [Equation 55] 
Normalization 

InNormal = In Tap /R Norma i (55) 

15 [0215] [Equation 56] 
Clipping 



In = < 



CLIP MAX In Nonnal > CLIP MAX 



CLIP_MIN In Nomial < CLEP_MIN (56) 
I n No m ai else 



[0216] In the case of 16 bits, CLIP MAX=32767 and CLIP 
MIN=32768 stand. 
20 In the case of 24 bits, CLIP MAX=8388607 and CLIP 

MIN=8388608 stand. 
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In the case of signed 8 bits, CLIP MAX=127 and CLIP 
MIN=128 stand. 

In the case of* unsigned 8 bits, CLIP MAX=255 and 
CLIP MIN=0 stand. 
5 [0217] <Step 16> 

Next, linear interpolation is carried out (F206) . 

The linear interpolation of one sample is carried 
out as follows. 

[0218] Minus minimum value of PCM 24bit*/ 
10 #define SRC PCM 24BIT MINUS-8388608 

Plus maximum value of PCM 24bit*/ 
#define SRC PCM 24BIT PLUS 8388607 

Maximum value of Alpha, Beta coefficients of Linear 
Upl28*/ 

15 #define SRC LINEAR MAX COEF 65535 

Rounding of 16 bit of linear interpolation*/ 
#define SRC LINEAR NORMAL 327 68 

Shift amount of 16 bit of linear interpolation*/ 
#define SRC LINEAR SHIFT 16 
20 [0219] Perform 1 sample of linear interpolation 

alpha, /*16 bit input*/ 
IN A, /*24 bit input*/ 
IN B, /*24 bit input*/ 
OUT Y/*24 bit output*/ 
25 lldata=alpha*IN B+ (SRC LINEAR MAX COEF-alpha) */IN 

A; 
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/*Rounding*/ 

if (lldata>=0) lldata+=(SRC LINEAR NORMAL); 
else lldata-=(SRC LINEAR NORMAL); 

/♦Normalize*/ 
5 Out Ya=(lldata»SRC LINEAR SHIFT); 

Clipping*/ 

Y(n)=(Out Ya<SRC PCM 24BIT MINUS) ? SRC PCM 24BIT 

MINUS 

: ((Out Ya>SRC PCM 24BIT PLUS) ? SRC PCM 24BIT 

10 PLUS : 

Out Ya) ; 
[0220] <Step 17> 

Update the counter value and the center of sample 
(F208) . 

15 The MasterCount and the sample number Readsample to 

be read are updated as follows. 

Add Count Value to control Master Counter */ 
MasterCount+=CountAdd; 

Remainder of decimal point or less*/1000 times */ 
20 Coun 1 1 0 0 0+=Coun tAmar i ; 

if(Countl000 >=1000) { 
MasterCount+=l ; 
Countl000-=1000 ; 
} 

25 InputOf f set+= ( (MasterCount) »23) &0x00000007 ; 

MasterCount=MasterCount&0x007f f f f f ; 
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[0221] <Step 18> 

Here, a judgment of whether or not the end 
condition is satisfied is carried out (F207) . When the 
condition is satisfied, the processing is ended. When it 
5 is not satisfied, the routine returns to the processing 
of F208 and F209. 

[0222] Next, an explanation will be given of the 
mounting method and a specific example of the sampling 
rate converter of FIG. 14, FIG. 21, and FIG. 22. 

10 [0223] <Mounting Method of Sampling Rate Converter 
(FIG. 14, FIG. 21, and FIG. 22)> 

When realizing the sampling rate converter shown in 
FIG. 14, it is realized by using the poly-phase 
configuration shown in FIG. 21. In the method explained 

15 below, the required smallest computation is satisfactory 
by selecting the required input with respect to the 
output and the poly-phase filter. 

[0224] FIG. 24 is a flow chart for explaining the 
mounting method of the sampling rate converter according 

20 to the present embodiment. 

The specific processing becomes as follows. Note 
that, for convenience of the explanation, C language-like 
expressions are used, but does not matter if the mounting 
means is hardware, software, or a mixture of hardware and 

25 software. 

[0225] The processing of FIG. 24 is the same as the 

68 



S05P0088 

processing of FIG. 23 except the processing of F209 in 
the processing of FIG. 23 is carried out after the 
processing of F206. 

The contents of the processing in the steps are the 
5 same as the contents explained with reference to FIG. 23, 
so overlapping explanations are omitted here. 
[0226] <Seventh Embodiment> 

As a seventh embodiment, an explanation will be 
given of an audio apparatus employing each above sampling 
10 rate converter . 

[0227] FIG. 25 is a block diagram showing an example 
of the configuration of an audio apparatus employing a 
sampling rate converter according to the present 
invention . 

15 [0228] The audio apparatus 300 has an input terminal 
301, LR separation circuit (DSB) 302, sampling rate 
converter (SRC) 303, attenuator (ATT) 304, mute circuit 
(MUTE) 305, and output terminal 306. 

[0229] As the sampling rate converter (SRC) 303, 
20 sampling rate converters 100, 100A, 200, and 200A to 200E 
explained with reference to FIG. 5, FIG. 14, FIG. 17, FIG. 
18, FIG. 21, and FIG. 22 as the above-explained first to 
sixth embodiments are applied. 

[0230] In this case, as exceptions to the above- 
25 explained explanation, the following processings are 
carried out. 
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[0231] End point processing is carried out. Namely, 
when a first filter computation is carried out, there is 
a shortage of Tap/2 samples. 0 is supplemented for the 
Tap/2 samples. 

5 Further, holding processing of previous data is 

carried out. As the processing peculiar to audio, it is 
necessary to copy the input data for the next filter 
computation after applying a filter with respect to Width 
as follows. 

10 [0232] [Equation 57] 

Src(i)=Src(Width+i) ,i=0,~, Tap-1 (57) 
[0233] In this audio apparatus 300, for example a PCM 
digital signal of 4 8 kHz input from the input terminal 
301 is subjected to a predetermined separating at a 

15 separation circuit 302, then the sampling rate is 

converted at the sampling rate converter 303, and the 
result is output as for example a signal of 44.1 kHz. 

Then, attenuation processing is applied at the 
attenuator 304, and further the PCM digital audio signal 

20 is output from the output terminal 306 via the mute 
circuit 305. 

[0234] According to the audio apparatus 300, since it 
has a sampling rate converter in which chessboard 
distortion is avoided, any pre-filter can be considered, 
25 any frequency point can be passed, and the amount of 

processing can be suppressed to the required lowest limit, 
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there are the advantages that the noise tolerance is 
improved, the deviation of the linear gain can be avoided, 
and an improvement of the processing speed can be 
achieved. 
5 [0235] <Eighth Embodiment 

As an eighth embodiment, an explanation will be 
given of an image processing apparatus employing a 
sampling rate converter. 

[0236] In the conversion of resolution of an image 
10 utilizing sampling rate conversion for the application of 
the filter processing using the linear phase FIR filter, 
use is made of a multi-rate filter using an interpolator, 
a decimeter, and a linear phase FIR filter as element 
technologies . 

15 [0237] In this case, the end point processing is 

carried out. As processings of both ends, processings 
such that the zero is supplemented first, the reflected 
(Mirror) data is input second, and pixels of both ends 
are held third are carried out. 

20 Further, the initialization processing is carried 

out. Specifically, after performing the processing up to 
the end of the image, CoefCount=0 is initialized. 
[0238] According to the present image processing 
apparatus, since it has a sampling rate converter in 

25 which chessboard distortion is avoided, any pre-filter 

can be considered, any frequency point can be passed, and 
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it is possible to suppress the amount of processing to 
the required lowest limit, there are the advantages that 
the image does not exhibit a lattice state, deviation of 
the linear gain can be avoided, and an improvement of the 
processing speed can be achieved. 

INDUSTRIAL APPLICABILITY 
[0239] According to the present invention, an 
amplitude characteristic that passes any frequency can be 
obtained, and high precision conversion can be 
accomplished without regard as to the cutoff frequency, 
therefore the present invention can be applied to for 
example conversion of the sampling frequency of audio and 
conversion of resolution such as magnification or 
reduction of the frame of images. 
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